Limits on Majoron-emitting double-beta decays of 136 Xe in the KamLAND-Zen experiment 
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We present limits on Majoron-emitting neutrinoless double-beta decay modes based on an exposure of 
112.3 days with 125 kg of 136 Xe. In particular, a lower limit on the ordinary (spectral index n — 1) Majoron- 
emitting decay half-life of 136 Xe is obtained as T°jf > 2.6 x 10 24 yr at 90% C.L., a factor of five more 
stringent than previous limits. The corresponding upper limit on the effective Majoron-neutrino coupling, using 
a range of available nuclear matrix calculations, is ((fee) < (0.8 — 1.6) x 10 . This excludes a previously 
unconstrained region of parameter space and strongly limits the possible contribution of ordinary Majoron emis- 
sion modes to Ov/3/3 decay for neutrino masses in the inverted hierarchy scheme. 

PACS numbers: 23.40.-s, 21.10.Tg, 14.60.Pq, 27.60. +j 



The search for neutrinoless double-beta (Ov/3/3) decay is the 
best probe of the Majorana nature of the neutrino known at 
present. The observation of this process would immediately 
imply total lepton number violation and the equivalence of the 
neutrino and the antineutrino, irrespective of the mechanism 
by which the decay is mediated [l]]. Although most current 
experimental efforts focus on the detection of Ov/3/3 decay 
mediated by light Majorana neutrino exchange, many other 
mechanisms have been proposed. Some exotic models |0, 
predict decays proceeding through the emission of massless 
Nambu-Goldstone (NG) bosons, referred to as Majorons. Pre- 
cise measurements of the invisible decay width of the Z boson 
in LEP [4] showed that traditional Majoron models require se- 
vere fine-tuning @]. However, a number of additional models 
have been proposed which avoid such fine tuning, including 
modes in which the Majoron can carry leptonic charge, and 
need not be a NG boson yg], or in which OV/3/3 decay proceeds 
through the emission of two Majorons (0]. These models pre- 
dict different shapes for the spectrum of the summed energy 
of the two emitted /3's. In this report we analyze the spec- 
trum obtained from a 38.6kg-year exposure of 136 Xe with 
KamLAND-Zen |8] to derive new limits on several of these 
decay modes. 

Table J] summarizes ten Majoron models lf/1, l9l-[TTll. which 
can be divided into two categories: (I) lepton number violat- 
ing models, and (II) lepton number conserving models. The 
table indicates also whether the corresponding Qvf3f3 decay is 



accompanied by the emission of one or two Majorons: 

(A,Z) -> (A,Z + 2) + 2e- + X ° , (1) 
(A,Z) -> (A,Z + 2) + 2e- + 2 X ° . (2) 

The main distinguishing features of the models are listed in 
the third, fourth, and fifth columns: whether the Majoron is 
a NG boson or not, its leptonic charge (L), and the model's 
spectral index (n), respectively. The spectral index is de- 
fined from the phase space of the emitted particles, G ~ 
(Q/3/3 — K) n , where Qpp is the Q value of the (3(3 decay and 
K is the summed energy of the electrons. A spectral index of 
n = 1 is denoted as ordinary Majoron emission, while n / 1 
is referred to as non-ordinary. Experimental searches for f3(3 
decay mediated by Majoron emission have been performed by 
Heidelberg-Moscow for 76 Ge H; by NEMO-2 and NEMO- 
3 for 100 Mo, 116 Cd, 82 Se, 96 Zr, 150 Nd, and 130 Te [QJ GJ- 
[11 ; by ELEGANT V for 100 Mo JH; and by DAMA for 
136 Xe ll2lll . Each of these searches reported limits for ordi- 
nary Majoron emission, while several also reported limits for 
non-ordinary emission modes. 



The experimental investigation reported here is based on 
data collected with KamLAND-Zen between October 12, 
2011, and February 9, 2012, and includes the data used 
in Ref. |§] plus an additional 34.7 days, for a total live- 
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TABLE I: Different Majoron emission models lffll9l 4Llll . and the KamLAND-Zen limits for the corresponding Majoron-emitting Ov/3/3 decay 
half-lives (T1/2) and effective Majoron-neutrino coupling constants ({g ee }) for 136 Xe at 90% C.L. The model notation in the first column 
follows Refs. (3] and fI3l . The third, fourth, fifth, and sixth columns indicate whether the Majoron is a NG boson or not, its leptonic charge 
(L), the model's spectral index and the form of the nuclear matrix element, respectively. 



Model 


Decay mode 


NG boson 


Li 


71 


Matrix element 


Results from this measurement 

Zi/a(yr) (See) 


IB 




No 





1 


M F - M GT [13, 14] 


> 2.6 x 10 24 


< (0.8 - 1.6) x 10~ 5 


IC 




Yes 





1 


Mp - M GT [13, 14] 


> 2.6 x 10 24 


< (0.8 - 1.6) x 10~ 5 


YD 


U^PPX X 


iNO 


A 

u 


1 

j 




> 4.5 x 10 23 


< 0.68 


IE 




Yes 





3 


M F „2 - Mq Tl0 2 [9J 


> 4.5 x 10 23 


< 0.68 


I1B 




No 


-2 


1 


Af F - Mgt [13, 14J 


> 2.6 x 10 24 


< (0.8 - 1.6) x 10~ 5 


IIC 




Yes 


-2 


3 


Mcfl [9J 


> 4.5 x 10 23 


< 0.013 


IID 




No 


-1 


3 


M F „2 - Mq Tuj 2 [9J 


> 4.5 x 10 23 


< 0.68 


HE 


0vj3p X °X° 


Yes 


-1 


7 


AW - M GT „ 2 [9J 


> 1.1 x 10 22 


< 1.2 


IIP 


QvPPx° 


Gauge boson 


-2 


3 


M CJi [9J 


> 4.5 x 10 23 


< 0.013 


"bulk" 




Bulk field 





2 




> 1.0 x 10 24 





time of 112.3 days. The target liquid scintillator (LS) con- 
tains (2.44 ± 0.01)% by weight of enriched xenon gas, ob- 
tained by augmenting the measurement by gas chromatog- 
raphy, (2.52 ± 0.07)% HI], with a much more precise esti- 
mate based on the total xenon weight and LS volume intro- 
duced into the detector during LS filling, (2.44 ± 0.01)%. Re- 
constructed energies and positions of muon-induced neutron- 
capture 7 events confirm that the detector status and data qual- 
ity were stable over the data-set. The fiducial volume for the 
reconstructed event vertices is defined as a spherical shape 
1.2 m in radius at the detector center, containing (125 ± 7) kg 
of 136 Xe. Other event selection criteria are applied as de- 
scribed in Ref. (Ht]. The systematic uncertainties are summa- 
rized in Table HTJ the dominant contribution comes from the 
fiducial volume uncertainty. The total systematic uncertainty 
on the /3f3 decay half-life measurement is 5.2%, slightly im- 
proved from the 5.9% of Ref. |8|] owing to the supplemental 
xenon concentration measurement. 

Figure Q] shows predicted energy spectra of 136 Xe decay 
for different values of the spectral index n, corresponding to 

2z//3/3 (n = 5), 0^x° (» = I, 2 , and 3), and 0^/3/3x°X° 
(n = 3 and 7). The spectra have been convolved with the 
KamLAND-Zen detector response function, including the en- 
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FIG. 1 : KamLAND-Zen visible energy spectra for different Xe 
decay modes, characterized by the spectral index n. The resolution- 
limited line [a = (6.6 ± 0.3)%/ ^/EXMeV)] at the Q value indicates 
the 0^/3/3 decay peak without Majoron emission. 



TABLE II: Estimated systematic uncertainties used for the 136 Xe f3/3 
decay half-life measurement. The overall uncertainty is 5.2%. 



Source 


Systematic Uncertainty (%) 


Fiducial volume 


5.2 


Enrichment factor of 136 Xe 


0.05 


Xenon concentration 


0.34 


Detector energy scale 


0.3 


Xe-LS edge effect 


0.06 


Detection efficiency 


0.2 


Total 


5.2 



ergy resolution and energy scale non-linearities. In this anal- 
ysis, the dominant contribution from 2^/3/3 decay is fit simul- 
taneously with a possible contribution from one of the Oi//3/3 
modes to obtain an upper limit on its rate. 

The energy spectrum of selected candidate events between 
0.5 and 4.8 MeV is shown in Fig. [2] The (3f3 decay rates are es- 
timated from a likelihood fit to the binned energy spectrum. In 
the fit, background contributions from external sources, from 
the 222 Rn- 210 Pb and 228 Th .208p b chainSi and from n C and 

10 C (muon spallation products), as well as the parameters of 
the detector energy response model are allowed to vary but 
are constrained by their independent measurement |18J. Back- 
grounds without independent measurements, namely 85 Kr, 
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FIG. 2: Energy spectrum of selected /3/3 decay candidates (data 
points) together with the best-fit backgrounds (gray dashed line) and 
2v/3{3 decay (purple solid line), and the 90% C.L. upper limit for 
Ovfifi decay and Majoron-emitting Oi/p/3 decays for each spectral in- 
dex. The red line depicts the sum of the 2^/3/3 decay and background 
spectra. Numerical results are reported in Table [[] The best-fit has a 
X 2 /d.o.f. = 100.4/87 for the full fit range 0.5 < E < 4.8 MeV. 



40 K, non-equilibrium 210 Bi, and the 238 TJ- 222 Rn and 232 Th- 
228 Th decay chains, are left unconstrained. As in Ref. JU, 
the background in the OV/3/3 region of interest, 2.2 < E < 
3.0 MeV, is fit to a combination of 110 Ag m , 88 Y, 208 Bi, and 
60 Co, constrained by the observed time variation of the event 
rate in that energy range, shown in Fig. [3] 

The additional exposure in this analysis yields an improved 
measurement of the 2^/3/3 decay half-life of 136 Xe. Setting the 
contributions from all Majoron-emitting modes to zero gives 
Tfj 2 = 2.30 ± 0.02(stat) ± 0.12(syst) x 10 21 yr, which is 
consistent with the previous result fH. The half-life limit for 
0^/3/3 decay also improves slightly to T®Y 2 > 6.2 x 10 24 yr at 
90% C.L. 

The event rate in the Ov/3/3 region of interest shows 
little time variation, limiting the allowed contribution of 
88y (T1/2 = 107 days) assuming its parent 88 Zr is ab- 
sent. There is little statistical power to distinguish 110 Ag" 1 
(T 1/2 = 250 days) and 208 Bi (T 1/2 = 3.7 x 10 5 yr). How- 



ever, the energy spectrum without ii0 Ag m is rejected by a 
X 2 test at more than 3er C.L., indicating a preference for a 
dominant contribution from this contaminant in the 0^/3/3 re- 
gion. We find that including a non-zero contamination of 88 Zr 
changes only the relative contributions of 88 Y and 208 Bi, but 
there is no impact on the other spectral components. 

The 90% C.L. upper limits for the different Majoron- 
emitting decay mode spectra are drawn in Fig. [2] and the 
corresponding half-life limits are listed in Table J] In partic- 
ular, the lower limit on the ordinary (spectral index n = 1) 



Majoron-emitting decay half-life is Tfj* > 2.6 x 10" yr at 
90% C.L., which is a factor of five more stringent than previ- 
ous limits ll2lll . Owing to the larger background at lower en- 
ergy, the sensitivity to higher spectral index decays is weaker. 
The limits on single- or double-Majoron emission can be 
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FIG. 3: Event rate variation in the energy regions (a) 1.2 < E < 
2.0 MeV (2vfiP window) and (b) 2.2 < E < 3.0 MeV (0^/3/3 
window). The fitted curves correspond to the expected variations 
for the hypotheses that all the events in the Oz//3/3 window are solely 
from one of the background candidates, 110 Ag m (solid line), 208 Bi 
(dotted line), or unsupported 88 Y (dot-dashed line). 



translated into limits on the effective coupling constant of the 
Majoron to the neutrino, {g e e) JUL using the relations 



1/2 



|( ffee )| 2 |M| 2 G forO^x , (3) 
T/2 = l(.9e e )| 4 |A/| 2 G forOw3/3x°x°. (4) 



The nuclear matrix elements M and the phase space factors 
G for the ordinary Majoron-emitting decay (n = 1) are taken 
from lfl3l[l4ll and $22% . respectively, while those for other de- 
cays are taken from H. From the half-life limit of the ordi- 
nary Majoron-emitting decay mode we obtain an upper limit 
of (g ee ) < (0.8 - 1.6) x 10~ 5 at 90% C.L., where the range 
of the upper limit corresponds to the theoretical range of the 
nuclear matrix elements lfl3l Q. This is the most stringent 
limit on (g ee ) to date among all /3/3 decay nuclei 1151 ll2L 
I2TI1 . The previous best limit from a laboratory experiment was 
fromNEMO-3 for 100 Mo: (g ee ) < (3.5 - 8.5) x 10~ 5 (H. 
Our new limit corresponds to more than a factor of 2.2 im- 
provement over this previous result. The limits on the effec- 
tive Majoron-neutrino coupling constant for 136 Xe for all in- 
vestigated Majoron-emitting 0^/3/3 decays are summarized in 
TableU 

Other limits on (g ee ) are available from geochemistry and 
astrophysics. Half-life limits on 128 Te from geochemical ex- 
periments can be interpreted as an effective coupling limit 
of (gee) < 3 x 10~ 5 lO, although the half-life determina- 
tions have been criticized l24ll and may require a down- 
ward correction by almost a factor of 3. The observation 
of neutrinos from SN1987A and of their time distribution 
indicates that Majoron emission does not play a dominant 
role in core collapse processes, allowing one to exclude the 
range 4 x 10~ 7 < (g ee ) < 2 x 10" 5 J^HH for the ordi- 
nary Majoron-emitting decay mode. While previous limits 
combined with the supernova data still allowed a gap re- 
gion of 2 x 10~ 5 < (g ee ) < 9 x 10~ 5 [il9], our new result 
completely excludes this region. The SN1987A limit there- 
fore significantly extends the KamLAND-Zen limit down to 
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(9ee) < 4 x 10~ 7 . Multiplying by the square-root of the ra- 
tio of phase space factors, one finds that this limit excludes 
the possibility that ordinary Majoron-emitting decay modes 
play a dominant role to light Majorana neutrino exchange for 
{ m pp) > 20meV. This range covers almost the entire (mpp) 
parameter space in the case of the inverted neutrino mass hi- 
erarchy JHf. 

In summary, we have reported new limits from KamLAND- 
Zen on Majoron-emitting 0^/3/3 decay half-lives for 136 Xe. In 
particular, for ordinary Majoron-emitting decay (spec- 

tral index n = 1), we obtained an improved limit on the effec- 
tive Majoron-neutrino coupling constant (g ee ) by more than 
a factor of 2.2 over earlier searches. Combined with limits 



on {g ee ) from SN1987A, the KamLAND-Zen result strongly 
disfavors a dominant contribution from ordinary Majoron- 
emitting decay modes for neutrino masses in the inverted hi- 
erarchy scheme. 
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